Interference of diffusing photons and level crossing spectroscopy 
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We show that a new interference effect appears in the intensity fluctuations of photons multiply 
scattered by an atomic gas of large optical depth b. This interference occurs only for scattering atoms 
that are Zeeman degenerate and it leads to a deviation from the Rayleigh law. The fluctuations 
measured by their variance, display a resonance peak as a function of an applied magnetic field. 
The resonance width is proportional to the small factor 1/6. We derive closed analytic expressions 
for all these physical quantities which are directly accessible experimentally. 
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We consider in this letter coherent multiple scatter- 
ing of photons propagating in cold atomic gases. For 
Zeeman degenerate atoms, we show that there is a new 
interference effect which affects the fluctuations of inten- 
sity (speckle pattern) of photons transmitted through a 
gas of large optical depth, but leaves unchanged the av- 
erage transmitted intensity. This interference shows up 
in a significant deviation from the Rayleigh law which 
states that for classical scatterers without internal struc- 
ture, the variance ST 2 = T 2 — T of the transmission 
coefficient T is simply related to its average value T by 
ST 2 = T . The averaging over configurations denoted 
by 7TT will be defined later. The interference involves the 
ground state Zeeman sublevels of the atoms. The depen- 
dence of this interference on magnetic quantum numbers 
suggests that it is sensitive to an applied magnetic field 
H . Our purpose in this letter is to show that indeed the 
amplified variance presents a resonance as a function of 
H around a crossing point, and the width of this reso- 
nance is 
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where I is the photon elastic mean free path through the 
atomic gas confined into a slab of width L, g is the Lande 
factor and [1q the Bohr magneton. Here, a is a constant 
of order unity to be determined later, that depends on 
the details of the atomic structure. All these features are 
obtained in the limit of diffusing photons, i.e., in a regime 
where L 3> I. Thus the narrowing of the resonance is in 
principle not limited, which might prove useful in level- 
crossing spectroscopic measurements. We shall see that 
the effect we present here shares some kind of analogy 
with the well-known Franken or Hanle effects [1, 2]. Nev- 
ertheless, we emphasize that both the underlying physical 
mechanisms and the quantities that are being measured 
are very different from these two effects [3, 4]. 

We consider the setup of Fig.l. A photon of polariza- 
tion e a is incident along a direction s a onto the atomic 
gas. It is detected in transmission with polarization if,, 
along Sf, after being multiply scattered. A time r later, 
a second identical photon is detected. We assume that r 



is short enough so that the atoms stay at rest between 
the two events. The same measurement is repeated after 
a time T 3> t, during which the scatterers move. The 
averaging over spatial disorder results from this motion. 
The transmitted intensity T is proportional to the prob- 
ability of a photon incoming along s a , to emerge along 
Sf,. Atoms are modeled as degenerate two-level systems. 
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FIG. 1: Experimental setup. A CW laser beam is incident 
onto atoms confined in a slab of width L. A detector placed 
along the direction Sb, records the intensity. During the time 
interval t, assumed to be short, the scatterers stay at rest. 
However, the two photons 1 and 2 experience different atomic 
internal configurations due to all other photons between them. 
The measurement is repeated after a time T>r. 

We denote by \m g ) the ground state Zeeman sublevels 
with total angular momentum j g , \m e ) are the excited 
states sublevels with total angular momentum j e , and m 
is the projection on a quantization axis. 

The average transmission coefficient T is obtained by 
squaring the sum of the scattering amplitudes, An R ' m ^ , 
corresponding to a given configuration {R, to}. Here {R} 
accounts for the spatial positions of all atoms, {to} is 
a notation for their internal Zeeman states both before 
and after scattering and the index n denotes one possible 
multiple scattering path. Then, 
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where • • • denotes a configuration average over both {R} 
and {to}. When averaging over {R}, all cross terms 
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n ^ n' vanish because of large fluctuating phase shifts, 

so that T ~ J2 n l^-rT^l 2 - This expression is the lead- 
ing approximation in the weak disorder limit fco^ S> 1, 
where ko is the photon wave number [5]. The two pho- 
tons detected at t = and t — t are separated by many 
undetected photons (sec Fig.l) which may change the 
internal states of atoms. Therefore, if {to} and {to'} re- 
fer to the atomic internal configuration seen by the two 
detected photons, then we can assume that there is no 
correlation between {to} and {to'}. 

Similarly, the correlation function of the transmission 
coefficients T and T' of the two detected photons is 
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As before, the averaging over {R} leaves only pairs of 
amplitudes having exactly opposite phase shifts. Thus, 
to leading order in weak disorder, the only non vanishing 
contributions involve two possible pairings of amplitudes, 
either i = j, k = I, which gives TT', or i = I, j = k, so 
that defining C 2 = TT 7 - TT, we obtain 
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This correlation function appears as products of two am- 
plitudes, that correspond to different internal configura- 
tions {to} and {to'}, but to the same scattering path 
i (or j). Most of multiple scattering paths i and j do 
not share common scatterers so that we can average 
A\ m ^ A\ m and A^ m * A^ m ^* separately, since these av- 
erages are taken upon different atoms, and finally, 
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Since generally {m} ^ {to'}, the interference occurs be- 
tween distinct Zeeman sublevels of the ground state, un- 
like the Franken or Hanle effects, where the interference 
involves distinct excited sublevels. This constitutes a new 
kind of interference, which originates from the fact that 
it is the correlation function rather than the average in- 
tensity, that is considered. 

In the theory of multiple scattering it is helpful to use a 
continuous description [5]. In this framework, one defines 
two Diffuson functions X>(* ,c ) by [6] 

T = j drdY'V {t) {ry ) and C= J drdr'V^ (r, r' 

(6) 

The two functions T>^ l,c ^ are obtained from an iteration 
equation (also called ladder diagram) whose structure is 
based on two elementary vertices y( l,c ), that describe the 
microscopic details of the scattering process. The itera- 
tion of the elementary vertices is written symbolically as 



where T>,V stand for T> < < % ' c \V^ l,c \ V accounts for a sin- 
gle scattering and VWV represents its iteration. The 
quantity W describes the propagation of the photon in- 
tensity between successive scattering events and it will 
be described later on. 

Generally, the elementary vertex is obtained by cou- 
pling two scattering amplitudes. It is given by 
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where the operator V(i',i) — £'* ■ d|TO e )(TO e |d • i results 
from the dipolar interaction energy — d.E between atoms 
and photons, d and E are respectively the atomic dipolc 
and electric field operators. The states \rm) are Zeeman 
sublevels of the atomic ground state, and \m e ), \m' e ) are 
those of the excited state. We have defined the energy 
difference HtOij = Ej — Ei and the photon frequency ui. 
We assume that the ground state Zeeman sublevels are 
equiprobable so that the corresponding density matrix 
reduces to the factor l/(2j g + 1). 

The elementary vertex that corresponds to the 

average intensity, is obtained by setting mi = TO3, to 2 = 
TO4, ii — £3 and £2 = £4 in (8). Up to a proportionality 
factor, is nothing but the differential cross section for 
this scattering process. Assuming a broad line excitation 
[8], we average over v, leading to [2] 
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where Bi 2 (m e ) — (to 2 |£2 ■ d|TO e }(m. e |d • EjJtoi). A mag- 
netic field removes the level degeneracy and leads to a 
Zeeman splitting, so that two kinds of terms appear in 
(9) depending on whether m e — m' e or m e 7^ m' e . Terms 
for which m e — m' e , are independent of magnetic field 
and give the incoherent scattering cross section. The 
terms m e ^ m' e depend on magnetic field and describe 
interferences between two distinct scattering amplitudes. 

For the vertex V^ c \ each one of the two coupled 
scattering amplitudes in (8) might belong to a distinct 
atomic configuration (see (5)), meaning that we must 
consider distinct couples of initial (|toi), |to 3 }) and final 
(lm.2), 1^4)) atomic states, as well as two initial (£1,^3) 
and final (£2, £4) polarization states. Summations over 
the quantum numbers to, result from averaging over ini- 
tial atomic states and from non detected final states. The 
vertex involves more interference terms than those 
already appearing in V^ . A non degenerate ground state 
leads immediately, using (8), to V (i) = V^ so that we 
recover the Rayleigh law, C 2 = TT [5]. Degenerate 
states produce additional interference terms in V^ so 
that C 2 > TT (see Fig. 2). An applied magnetic field 
removes the degeneracy and therefore affects the inter- 
ference pattern. 

We now calculate the Diffusons T> given by the iter- 
ation (7). We first decompose the various terms into 
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FIG. 2: Example of a term, which, for the transition j g — 
1 — » j e = 2, contributes to V' c ', but not to V^. The solid 
and dashed lines refer respectively to the two quantum ampli- 
tudes that appear in (8). For large enough magnetic field, this 
contribution vanishes since Zeeman splitting takes it far from 
resonance. 



components in a standard basis, 
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Likewise, Eq.(7) acquires a tensorial form, 

£>a/3 :7 <5 = Va/3,j8 + W ^ T> a 0^ v 'P^ t p Cr Vp (7tl S ■ (H) 



Here W describes the scalar part of W and V a p,-y6 — 
{(S aj — (— ) 7 s Q s_ 7 ) (Sps — (— ) l3 s-p§s)) accounts for the 
polarization dependent part. This follows at once by 
noticing that after being scattered by an atom, the two 
outgoing photon amplitudes propagate with a wavevec- 
tor s = k/fco, random in direction but identical for both, 
and with two different polarization components. Since s 
is random, the intensity propagation is averaged (•••) over 
photon wavevectors direction. The term 5^ — (— )' y s AI s_^ 
expresses transversality. The two tensors V a f3,j8 and 
V Q( 3 i7 5 can be written as 9 x 9 matrices. The iteration 
(7) rewrites V = (1 + WW + (WW) 2 + • • -)V. We 
now use the spectral decomposition theorem to expand 
W = J2 K UkT^ k \ where uk's are the eigenvalues of 
W and the T^'s define an orthonormal set of (gener- 
ally) 9 projectors [7]. Then, with the help of (11), we 
find 
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where A = (2j e + l)/3(2j g + 1) and q (with q = |q|) is 
the Fourier variable of the difference R = r' — r between 
the two endpoints of a multiple scattering sequence. The 
r.h.s in (13) is obtained by using the diffusion approx- 
imation (i.e. ql <C 1), so that W(q) ~ ^ (1 - q 2 l 2 /3), 
where D = cl/3 is the photon diffusion coefficient [5]. We 
identify the set of characteristic damping rates 
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The term |A is the total cross section conveniently nor- 
malized. According to the values of u^ c \ we identify 3 
kinds of modes. A positive 7^'^ describes an exponen- 
tially damped mode. A vanishing 7^'°' corresponds to an 
infinite range stable mode which ensures energy conserva- 
tion, and a negative damping rate describes an amplified 
mode. The largest eigenvalue of V^V is = 2A/3, 
thus leading to one stable mode. For degenerate scatter- 
ers (jgjje > 0), and without magnetic field, has one 
amplified mode, whose occurrence results from the fact 
that all terms that contribute to contribute also to 
[7]. However, there are interference terms that con- 
tribute to only (Fig. 2). Therefore, the largest eigen- 
value of V^V becomes greater than |A, thus making 
the corresponding damping rate negative. 
We now rewrite (12) in real space, 
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where Y^'^ are two angular functions that depend on the 
incoming and outgoing polarizations e a and £& [7]. The 
scalar Diffuson propagator D(r, r',t) obeys a diffusion 
equation whose solution for a slab geometry is well known 
[5] and leads for (6) to 
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where we have defined L 
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transmission coefficient T is given by the same relation 
(16) provided we replace (L { ^,Y^ c) ) by (L^ ] ,Y^ } ) de- 
fined accordingly 

The dominant contribution to the average intensity T 

(i) 

is the stable, energy conserving mode Jo = 0- The two 
other modes have positive damping rates and are negligi- 
ble compared to this stable mode. They express photon 
depolarization in multiple scattering. The stable mode 
leads to T oc l/L = 1/naL. Here n is the density of scat- 
terers and a is the single scattering total cross section. 
However, the total cross section is independent of the 
magnetic field H [2]. This can be understood as follows. 
Starting from (9), the outgoing polarization dependent 
part is a sum of terms like X)e 2 ±k( TO2 |s2 ' d|m e ) (m' e \e2 ■ 
d|m 2 ) = ]Ck( £ 2a£2/3 + e2 a £' 2 p){m2\d a \m e )(m' e \d-p\m 2 }, 
where k is the outgoing wave vector, and (a, [3) are the 
standard components of d that contribute to the transi- 
tions. Integrating over k imposes a — —(3, which implies 
m e = m' e . Thus the interference terms in (9) does not 
contribute to a, which is therefore independent of H . As 
a consequence, T is also independent of H . 

The intensity correlation function C is dominated by 
the amplified mode driven by the negative damping rate 
7q C \ The integral in (15) is cutoff by t max = C/c where 
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C = cL 2 /D is the longest path of a diffusing photon. 
Eq.(16) thus leads to [7] 
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where X = L/L^ and b = L/l is the optical depth. 
This expression is displayed in Fig. 3 as a function of the 
dimcnsionless magnetic field s = g[ioH/frT. It retains the 
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FIG. 3: Plot of C given by (17) as a function of the dimen- 
sionless magnetic field s. We denote by A its FWHM. This 
plot corresponds to the j g — 1 — » j e = 2 transition and 6 = 5. 



interference to a magnetic field is significantly enhanced. 
This could be used towards more precise measurements 
in level crossing spectroscopy, in the limit of dense atomic 
gases where multiple scattering cannot be neglected any- 
more. 
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FIG. 4: FWHM A of C plotted versus 1/b. The points are 
obtained from (17) in which fi results from a direct nu- 
merical calculation of Jq C \s). The slope of this linear be- 
havior is in good agreement with the predicted expression 
2y^ln 2//2 ~ 0.76 obtained for the transition j g = 1 — » j e = 2. 



shape of a resonance curve whose FWHM is given by (1). 
To derive it, we expand V^ c \s) to first order in s, leading 
for the amplified mode to u^\s) ~ u^\o) — [3s 2 where 
[3 is a constant that depends on the specific scattering 
atom. According to (14), the corresponding damping 
rate becomes, 

7^ ) W-^ ) (0) + ^, 2 (18) 



By rewriting X = by 1 2 ^/q (s)/ D, and making use of 
(18), we obtain X = by/\fo — /2s 2 1, where the two con- 
stants /o and /2 are given by fo = (2A/u^) — 3 and 
fi = 2 A [3 /u^ 2 . For large enough optical depth b = L/l, 
i.e., in the diffusive regime where expression (17) ap- 
plies, the FWHM A behaves linearly with 1/6 as shown 
in Fig. 4. The corresponding slope is easily calculated 
from (17) and restoring units, we obtain for AH the ex- 
pression (1) with a = 2y4n2//2. 

To summarize, we have shown that in multiple scatter- 
ing, interference of diffusing photons scattered by atoms 
in the presence of a magnetic field near a level cross- 
ing (or close to zero field) shows up as a resonance peak 
in the intensity correlation function C. Its width, given 
by Eq.(l), is inversely proportional to the optical depth 
b = L/l. The diffusive regime corresponds to large val- 
ues of b (typically b ~ 10 2 ), so that the sensitivity of the 
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